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C ' J Abstract. The main objective of this article is to study both dynamic and 

^SJ structural transitions of the Taylor-Couette flow, using the dynamic transition 

theory and geometric theory of incompressible flows developed recently by 
the authors. In particular we show that as the Taylor number crosses the 
critical number, the system undergoes either a continuous or a jump dynamic 
^^■ t transition, dictated by the sign of a computable, nondimensional parameter 

R. In addition, we show that the new transition states have the Taylor vortex 
type of flow structure, which is structurally stable. 
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1. Introduction 



The study of hydro-dynamic instability caused by the centrifugal forces originated 
from the famous experiments conducted by [13] in 1923, in which he observed and 
studied the stability of an incompressible viscous fluid between two rotating coaxial 
cylinders. In his experiments, Taylor investigated the case where the gap between 
the two cylinders is small in comparison with the mean radius, and the two cylinders 
rotate in the same direction. He found that when the Taylor number T is smaller 

i2; than a critical value T c > 0, called the critical Taylor number, the basic flow, called 

the Couette flow, is stable, and when the Taylor number crosses the critical value, 

£N) the Couette flow breaks out into a radially symmetric cellular pattern as in Figure 

in EU 

(~) There have been extensive studies for the Taylor problem from both the mathe- 

matical and physical point of view; see among many others, [16j [HI [IJ [2] • Over the 
years, the Taylor problem, together with the Rayleigh-Benard convection problem, 
^. has become one of the paradigms for studying nonequilibrium phase transitions and 

pattern formation in nonlinear sciences. 

The main objective of this article is to address the dynamic transition of the 
Taylor-Couette flow, and and study the formation and stability in its structure of 
the Taylor vortices. The main technical tools are the dynamical transition theory 
and the geometric theory for incompressible flows, both developed recently by the 
authors; see [SJ [TT] and the references therein. 

The main philosophy of the dynamic transition theory is to search for the full set 
of transition states, giving a complete characterization on stability and transition. 
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Figure 1.1. Couette Flow and 
Taylor vortices 




The set of transition states is represented by a local attractor. Following this phi- 
losophy, the dynamic transition theory is developed to identify the transition states 
and to classify them both dynamically and physically. One important ingredient 
of this theory is the introduction of a dynamic classification scheme of phase tran- 
sitions. With this classification scheme, phase transitions are classified into three 
types: continuous (Type-I), jump (Type-II) and mixed (Type-Ill). The dynamic 
transition theory is recently developed by the authors to identify the transition 
states and to classify them both dynamically and physically; see above references 
for details. The theory is motivated by phase transition problems in nonlinear 
sciences. Namely, the mathematical theory is developed under close links to the 
physics, and in return the theory is applied to the physical problems, although more 
applications are yet to be explored. With this theory, many long standing phase 
transition problems are either solved or become more accessible, providing new 
insights to both theoretical and experimental studies for the underlying physical 
problems. 

For simplicity, we focus in this article on the z-periodic boundary condition, 
which is an approximate description for the case where the ratio Lj(ri — r\) between 
the height L and the gap r 2 — r% is sufficiently large. We remark that similar 
results hold true as well for other type of boundary conditions, as well as for three 
dimensional perturbations (in the narrow-gap case); we refer the interested readers 
to [11] for further details. 

The main results obtained are as follows. 

First, we show that the system always undergoes a dynamic transition as the 
Taylor number T crosses the critical Taylor number T c . The types of the transition 
can be either continuous (Typc-I) or jump (Type-II), and are dictated precisely by 
the sign of a nondimensional parameter R, given completely by the first eigenvec- 
tors, the ratio of the angular velocity of the outer and inner cylinders /x, and and 
the ratio of the radii of the inner and outer cylinders rj. 

Second, when R < 0, the transition is continuous, and the critical exponent 
of the phase transition, i.e., the exponent in the expression of bifurcated solutions, 
is f3 — 1/2. Moreover, there is only one critical Taylor number T c such that the 
secondary flow tends to the basic flow (Couette flow) as T — >• T c . 
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Also, for the narrow-gap case, the parameter R defined by (3.231 is negative: 
R < 0, provided the two coaxial cylinders rotating in the same direction, including 
the case where the outer cylinder does not rotate. 

Third, when R > 0, the transition is a jump transition, leading to more drastic 
changes, coexistence of metastable states, and potentially more chaotic/turbulent 
behavior. In particular, there are two critical Taylor numbers T c and T* with 
T* < T c . When T* < T < T c , the system has two metastable states E , the 
trivial Couette flow and S T , a local attractor away from the Couette flow. When 
T > T c , the solution always moves away from the basic Couette flow to a more 
chaotic/turbulent regime, represented by the local attractor E T . 

Fourth, the theoretic analysis carried out in this article shows that a street of 
vortices appear in the secondary flow for the narrow-gap case with /i — > 1. Thus 
the theoretic results are in agreement with the Taylor experiments. 

The article is organized as follows. The partial differential equation model and 
the set-up are given in Section 2, and the main dynamic transition theorems are 
given in Section 3. Explicit expressions of the parameter R for determining the 
types of transitions are further discussed in Section 4. The formation and structural 
stability of the Taylor vortices are addressed further in Section 5, and the main 
theorems are proved in Section 6. 

2. The Taylor Problem 

2.1. Couette flow and Taylor vortices. Consider an incompressible viscous 
fluid between two coaxial cylinders. Let r\ and r 2 (r 2 > r{) be the radii of the 
two cylinders, Q\ and f2 2 the angular velocities of the inner and the outer cylinders 
respectively, and 

(2.1) n = n 2 /n x , r i = r 1 /r 2 . 
The nondimensional Taylor number is defined by 

(2.2) T=^^, 



where v > is the kinematic viscosity, and h is the vertical length scale. 

There exists a basic steady state flow, called the Couette flow. In the cylindrical 
polar coordinate (r, 9, z), the Couette flow is defined by 

, {u r ,u ,u z ,p) = f 0,V(r),0,p / -V 2 (r)dr 

V(r) = ar + b/r, 

where (u r ,ug,u z ) is the velocity field, p is the pressure, and a, b are constants. It 
follows from the boundary conditions that 

V( ri ) = n iri , V{r 2 ) = Q 2 r 2 , 

and the constants a and b in (|2.3|) are given by 



a = -iliTj— =-, fr = S2i— =-. 

1 — r] z 1 — i] z 



where /i and 77 are given by (2.1 ) 
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Based on the Rayleigh criterion, when p > rj 2 , the Couette flow is always stable 
at a distribution of angular velocities 

il(r) = a + b/r for t\ < r < r-2- 

However, when p < rj 2 , the situation is different. As in the Taylor experiments, 
consider the case where the gap r 2 — T\ is much smaller than the mean radius 

r = 1/2(7-! + r 2 ), namely, 

r 2 -ri < (ri+r 2 )/2, 



and the two cylinders rotate in the same direction. If the Taylor number T in (2.2) 



satisfies T < T c , then the Couette flow (2.3) is stable, and if T c < T < T c + e for 



some s > 0, a street of vortices along the 2-axis, called the Taylor vortices, emerge 
abruptly from the basic flow, as shown in Figure |1.1| and the corresponding flow 
pattern is radically symmetric and structurally stable. 

When the gap r 2 — r 2 is not small than tq — 1/2 (n + r 2 ), or when the cylin- 
ders rotate in the opposite directions, the phenomena one observes are much more 
complex; see p] for details. 

Hence, in this section we always assume the condition 

(2.4) rj 2 > p > 0. 

2.2. Governing equations. The hydrodynamic equations governing an incom- 
pressible viscous fluid between two coaxial cylinders are the Navier-Stokes equa- 
tions. In the cylindrical polar coordinates (r, 0, z), they are given by 



(2.5) 



du r u 2 ( 2 due u r \ 
¥ + (u . VK - 7 = ^A« r -^ ¥ --5-J- 
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where v is the kinematic viscosity, p is the density, u — (u r ,ue,u z ) is the velocity 
field, p is the pressure function, and 
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Then it is easy to see that the Couette flow (2.3) is a steady state solution of 



(2.5). In order to investigate its stability and transitions, we need to consider the 



perturbed state of (2.3) 



u r ,u e + V(r),u z ,p + p / -v 2 (r)dr. 
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The perturbed equations read 



du r 

Ik 



+ (u ■ V)u, 



_ n 

r 
2V(r) 



2 du 6 



I dp 
p dr 



r 



-Ug 



V(r) du r 
r ~W 



(2.6) 



2 di 



dug , , _s . "eur , ,, 



no 



rp cw r r at* 






+ (u ■ V)u z 



vAu z 



d(ru z ) d(ru r ) due 



dz 
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do 



I dp V du z 
p dz r d9 

0. 



To derive the nondimensional form of equations (2.6 1, let 

(x, t) = (hx\ h 2 t'/v) (x = (r, rO, z)), 

(u,p) = (vu'/h, pv 2 p /h 2 ) (u = (u r ,ug,u z )). 



Omitting the primes, we obtain the nondimensional form of (2.6) as follows: 

..2 
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where T is the Taylor number as defined in (2.2). 
The nondimensional domain for (2.7) is 

n=(Ji,fa)x(0 > 27r)x(0 J i), 

where l- L = ri/h (i — 1, 2), and L is the height of the fluid between the two cylinders. 
The initial value condition for (2.7) is given by 

(2.8) u(r,$,z,0)=u (r,e,z). 

There are different physically sound boundary conditions. In the 0-direction it is 
periodic 

(2.9) u(r,0 + 2kn,z) =u{r,0,z), Vfc E Z. 
In the radical direction, there is the rigid boundary condition 

(2.10) u = {u z ,u r , ug) — 0, at r — li,l 2 - 
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At the top and bottom in the ^-direction [z = 0,L), either the free boundary 
condition or the rigid boundary condition or the periodic boundary condition can 
be used: 

Dirichlet Boundary Condition: 

(2.11) u = (u r , ug, u z ) = at z — 0, L; 

Free-Slip Boundary Condition: 

du r dug 

dz dz 
Free-Rigid Boundary Condition: 

- n ^^ — ^i 
(2.13) Uz ~ '~dz~~~dz~ z 

u = (u z ,u r ,u g ) = 
Periodic Boundary Condition: 



(2.12) 



0, 



at z = 0, L; 

— at z = L, 
at z = 0; 



(2.14) 



u(r,8 7 z + 2kL)=u(r,9z) Vfc G Z. 



3. Dynamic Transitions 



3.1. Functional setting. We now study the Taylor problem (2.7) with the z- 
periodic boundary condition (2.14) and with axisymmetric perturbations. Assum- 
ing that the equations ( 2.7 1 are independent of 9, and taking the length scale h = r% 
in the nondimensional form, we obtain 



du z dp 

— — = Au z - (u- V)u z , 

dt dz 



(3.1) 
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where A = \T, T is the Taylor number, and 
T= 4rf^(l-^)V ^ 

V 2 {\ — Jlj 2 ) 2 

1 -fX 
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(u- V) =u r — +u z — . 
or oz 

The nondimensional domain is M — (r), 1) x (0, L), and the boundary conditions 



v 2 - 


2 / 2 


A = 


a 2 19 a 2 

<9r 2 r dr dz 2 



take ( |2.10| and \2.U\, i.e., 

(3.2) 



u = (v, z ,u r ,ue) = at r = r/, 1, 

u is periodic with period L in the z — direction. 

The initial value condition is 

(3.3) u = Uo(r, z) at t = 0. 
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For the Taylor problem (3.1)-(3.3), we set 

_ J _ ,_ v T 2{i\, f \3 div(ru) = 0,u r = at r = 77, 1, and 
1 ^ ' ' ^ ' u is L — periodic in the z — direction 

Hi = U e H 2 (M) 3 f]H | u satisfies (3.2 )| , 

and the inner product of H is defined by 



{u,v)h — / u ■ vrdzdr. 

Let the linear operator L\ = —A + XB : Hi — > H and nonlinear operator 
G : Hi -)• P be defined by 



Aw 



-P ( Aw 2 , I A - — )u r , ( A - — I u g 



1 



— - K | U S ,KU r 



(3.4) Pit = P ( 

G(u) = -P( (u'VK(S'V)a 



— , (U- V)Ug H 

r r 



where P : L 2 (M) 3 — > H is the Leray projection. Thus the Taylor problem (3.1) 
(3.3) is rewritten in the abstract form 

dn 

(3.5) d£ 

u(0) = u Q . 



L\u + G(u), 

-- u . 
For simplicity, let G : Hi — >• H be the corresponding bilinear operator defined 



by 



G(u, v) = -P ((u ■ V)v z . (u ■ V)u r -—,(«• V)v e + ^ 
V r r 



Then it is easy to see that 

(3.6) (G(u, v),w) H = ~(G(u, w),v) H . 

3.2. Eigenvalue problem. To study the phase transition of the Taylor problem 
(3.1 )-(3.3) it is necessary to consider the eigenvalue problem of its linearized equa- 
tion. The associated eigenvalue equation of (3.5) is as follows: 

(3.7) L x u = -Au + XBu = (3(X)u, 



and the conjugate equation of (3.7) is given by 
(3.8) 



Liu* 



A*u* +XB*u* =/3(A)u* 



The equations corresponding to (3.7) are as follows 

dp 



Au z -^f= /3(X)u z 
oz 



(3.9) 



dp 



A ^)u r + x( — -k\u 9 -^ = /3(A)u r , 



1 



Ug + Xnu r = (3(X)ug, 



div(ru) = 0. 
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The equations corresponding to ( |3.8[ ) are given by 

dp* 
dz 



(3.10) 



A 
A 



/3(A)< 



1 



u* + Xnut 



dp* 
dr 



u* e + X 



/?(AK, 



(3(X)u* g 



div(rw*) = 0. 



Both (3.9) and (3.10) are supplemented with the boundary condition (3.2). 

We start with the principle of exchange of stability (PES). It is known that 
for each given period L, there is a Aq = \q(L) such that the eigenvalues /3j(X) 
(j = 1, 2, • • • ) of ([3T9]) with ([3T2]) near X = X* Q satisfy that /3i(A), • • • , fi m {X) (m > 1) 
are real, and 



(3.11) 



A (A) 



< if A < Aq, 
= if A = A^ 



for 1 < i < to, 



Re/^(Ao)<0 fori>m + l. 

In addition, there is a period L' > such that 
(3.12) A = Xo(L') = mmAo(L). 

L>0 



Thanks to [THIHI], for /i = SI1/O2 > the multiplicity m = 2 in (3.11) at Ao = 
Ao(L'); see also (3E]. 

In this section, we always take L' as the period given by (3.12), and define the 
following number as the critical Taylor number: 

Tc = Xl(L'). 

For simplicity, omitting the prime we denote L' by L. 

By (3.11) and (3.121, to verify the PES it suffices to prove that for A > Ao, 

(3.13) ft(A) > VI < i < m. 



To this end, we need to derive the eigenvectors of (3.9) and (3.10) at /3,(Ao) = 
(t = l,2). 

It is readily to check that the eigenvectors of (3.9) with (3.2) corresponding to 
A(Ao) = (i = 1, 2) are given by 

(3.14) ipi = (ip z ,ip r ,ip0) = (— sina^; D*h(r), acosaz h(r),cosaz <p(r)), 

(3.15) ipi = (ipzjipnipe) = (cosaz D t h(r),asmaz h(r), sinaz f{v)), 
where (h(r),ip(r)) satisfies 



(DD* - a 2 fh = a 2 X Q [ -=• — «) y? 3 



(3.16) 



and 



a 2 )<p 



{DD* 
(h,Dh,ip) =0 



-Aqk/i, 



D = 



d 
dr ' 



D* = 



at r = 77, 1, 

2- 



dr 
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The dual eigenvectors of (3.101 with (3.2 1 read 

(3.17) ip\ = (ip*, ip*, ip*,) = (— sinaz D*h* (r), a cos a z h*(r),cosaz ip*(r)), 

(3.18) ip\ = {ip* , ip* , ipl ) = (cosaz D« h* (r) , a sin az h*(r), sinaz f*(r)), 
where (h*,ip*) satisfies that 

(DD* ~ 2 ^ 2 



) h* = XoKtp* , 
(3.19) (DD, - a 2 )ip* = -a 2 \ f \ - k ) h 

(h*,Dh*, l f*) = 



at r = rj, 1. 



The following lemma shows that the PES is valid for the Taylor problem (3.1) 



(3.3) with [i > 0. 



Lemma 3.1. If [i > 0, then the first eigenvalues Pi{\) (1 < i < m) of {3.9) are real 
with multiplicity m — 2 near A = Ao = \Tc, and the first eigenvectors at X = Xq are 
given by {3.14) and (3.15). Moreover, the eigenvalues (3j(X) (j = 1,2,---) satisfy 
the conditions (5.4) and (5.5) at X = A , i.e., the PES holds true at the critical 
Taylor number T c . 



Proof. We only need to prove (3.13). By Theorem 2.1 in [jj], it suffices to verify 
that 



(3.20) 



(BVi.ViWO, {B^l) H ^Q. 



We infer from ((3_4j), (|3_14f, (j3_15j), fl3.17[ ), and ( |3.18[ ) that 
(3.21) (Bi'^i^H 

1 



L pi 
/ r 

2 Jr, 



K Ipelp* + K1p r tpg 



dzdr 



1 



— k I hip* + Kiph* 



dr. 



Since /i > 0, by ( 2.4 ), we have < k < 1 and -3 — k > for rj < r < 1. On the other 
hand, we know that the first eigenvectors (h(r),ip(r)) of (3.16) and (h* (r) , ip* (r)) 



of (3.19) at A = A are positive; see [15l l3"l IT4"] : 
(3.22) h(r)>0, (p(r)>0, h*(r)>0, 



ip*(r) > V77 < r < 1. 



Thus (3.20) follows from (3.21) and (3.22). The proof is complete. 



□ 



3.3. Phase transition theorems. Here we always assume that the first eigenvalue 



of (3.9) with (3.2) is real with multiplicity m = 2, i.e., the first eigenvalue Ao of 



(3.16) is simple, and the PES holds true. By Lemma 3.1 this assumption is valid 



for all /i > and < 77 < 1. 

Let ipx and ip* be given by (3.14) and (3.17). We define a number R by 
1 



(3.23) 



R 



(iPi,iI>*)h 
where $ € Hi is defined by 

(3.24) (A 



(G($>i),V'i)h + (G(Vi,$),V'i)h] ! 



A B)$ = G(^i,Vi 



Here the operator A, B and G are as in (3.4). The solution $ of (3.24) exists 



because G(i/),ip{) is orthogonal with ip* and ip* in H. 
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The following results characterize the dynamical properties of phase transitions 
for the Taylor problem with the z-periodic boundary condition. 



Theorem 3.1. If the number R < in \3.23), then the Taylor problem [3.1) 



(3.3) has a Tyep-I (continuous) transition at the critical Taylor number T — T c or 
X = \q, and the following assertions holds true: 

(1) When the Taylor number T <T C or A < Ao, the steady state u = is locally 
asymptotically stable. 

(2) The problem bifurcates from (u, A) = (0, Ao) (or from (u,T) — (0, T c )) to 
an attractor A\ homeomorphic to a circle S 1 on Ao < A, which consists of 
steady states of this problem. 

(3) Any u £ A\ can be expressed as 

u=\(3 1 (X)/R\ l ^v + o(\(J 1 \ 1 ^), 
v — xipx + yipx, 
x 2 + y 2 = l, 



where ipi,ipi are given in {3.14) and (3.11). 

(4) There is an open set U C H with £ U such that A\ attracts U\T, where 
r is the stable manifold of u = with codimension two in H. 

(5) When 1 — jj, > is small, for any u £ U \ (T\JH), there exists a time 
to > such that for any t > to, the vector field u(t,tio) = (u z ,u r ) is 
topologically equivalent to one of the patterns shown in Figure \3.1\ where 



u = (u(t,uo),ug(t,uo)) is the solution of \3.1)-(3.3), and 



L r l 



H = {u = (u z ,u r ,ug) £ H\ / ru z drdz = 0}. 

JO J-q 

(6) When 1 — \i > is small, for any uq £ (Uf]H) \ T, there exists a time 
to > such that for any t > to^u^t, Uq) = (u z ,u r ) is topologically equivalent 
to the structure as shown in Figure \S~^ 

Theorem 3.2. For the case where R > 0, the transition of the Taylor problem 



(3.1)-(3.3) at T — T c is of Type-II. Moreover, the Taylor problem has a singularity 

separation at T* < T c (A* < Ao). More precisely we have the following assertions: 

(1) There exists a number X* (0 < A* < Ao) such that the problem generate 

a circle E* = S 1 at X = X* consisting of singular points, and bifurcates 

from (£*, A*) on X* < X to at least two branches of circles E^ and E2, each 

consisting of steady states satisfying 

lim Ei - {0}, 



-.A 



dist ( T,n, 0) = mm \\u\\ H >0 at A = A 

uesj 



see Figure [X 

(2) For each X* < A < Ao, the space H can be decomposed into two open sets 
Ul and J7 2 A : H = tj£ + C/ 2 A with Uf f| C/ 2 A = 0, E A C dU£ f| dU$ such that 
the problem has two disjoint attractors A\ and _A A : 

A x x = {0} c U$, E A c A\ c (7 A , 

and Af attracts C/ A (i = 1,2). 
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(a) 



(b) 



Figure 3.1. Taylor 
vortices with a 
cross-channel 
flow. 







-^ r 



Figure 3.2. Taylor vortices without a 
cross-channel flow. 



(3) For Ao < A, the problem has an attractor A x satisfying 
lim A x = A Xo , dist(A x 7 0) > VA > A , 

A-s-Ao 

and A x attracts H \T\, where T\ is the stable manifold of u — with 
codimension m\ > 2 in H . 



12 
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Figure 3.3. Singularity 
separation of circles 
consisting of steady 
states at A = A*. 




4. Explicit expression of the parameter R 



4.1. General case. The parameter R defined by (3.23) and (3.24) can be explicitly 
expressed in the following integral formula: 

R = - 



oo 



rhip* — — dr 
dr 



J f r((4> ■ v)^ + (i> ■ V)4>A1 + (4> ■ V)A€ 

Jri V 



+ 



drdz 



where -0 = (i/) 2 ,ip r ),ij>i = {ip z ,?p r ,ipg),ipl = {ip* z , ip* , tp* e ) are given by ( |3.14[ > and 
§37} , 

(^i,-0i)h= / / r(ip z ip* z +ip r ip* +if)gipg)drdz, 

Jo J?) 

and 0o, = (<j> z ,<j) r ,<f>d) satisfy 
D£U 



a(ipD*h + hDip -\ — v 3 ' 1 ); 
r 



^0|r=77,l = 0, 

dz 

- ! A - -j J <?V - A ( -2 - « J <t>9 + ■£ = -l/2co&2azH 2 {r), 



A(f> z 



■l/2sm2azH 1 (r), 



- A- 



1 



<9r 



(j)g — XoK(j) r = — l/2cos2aziJ3(r), 



div(r0) = 0, 



Here iJj(r),i — 1,2,3, are as in (6.9). 



4.2. Narrow-gap case. We consider here the case where the gap r 2 — n is small 
compared to the mean radius r§ = (n + r-i)/2 with /x > and with axisymmetric 
perturbations. This case is the situation investigated by (T3] in 1923. 

We take the length scale h — r 2 — r± . Then the narrow gap condition is given by 



(4.1) 



1 =r 2 -r 1 < (n +r 2 )/2. 
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Under the assumption (4.1), we can neglect the terms containing r n (n > 1) in 

= Vt(i 



(2.7). In addition, by (4.1) we have 



H 1 - fi r\ 



T 

1 — rj 2 1 — rf 2 r 2 



I — H r — r\ 

1 — rj 2 r 2 



T(l-(l- M )(r-n)). 



Let 

(4.2) 



a 



2 

rr_-ji 
I-77 2 ' 



Replacing ug by y/aug, and assuming the perturbations are axi-symmetric and are 
independent of 9, we obtain from (2.7): 



du. 



+ (u ■ V)u z 



(4.3) 



dt 

du r 

Ik 

du e 

dt 

du T du. 

— - H 

dr dz 



dp 
dz' 

dp 



Au z 

+ (u ■ V)u r — Au r — ^ + VaT(l - (1 - u)(r - ri))u$, 
or 



where 



u ■ V)ug = Aug + VaTu r 
0. 



d 2 d 2 ,_ „, d d 
(tt .V) -«,- + «,_. 



9r 2 9z 2 

In this case, the spatial domain is M = (r 1; r\ + 1) X (0, L). For convenience, we 
consider here the Dirichlet boundary condition 

(4.4) u\ dM = 0. 

The initial value condition is axisymmetric, and given by 

(4.5) u = u n (r, z) at t = 0. 



The linearized equations of (4.3) read 

dp 



-At 



dz 

dp 



= 0. 



(4.6) 



where A 



- Au r + — = Xu e - A(l - ^)(r - ri)ug, 

or 

— Aug = \u r , 

du r du z 

— - H = 0, 

dr dz 

'aT, T is the Taylor number given by (|2.2[) 



Let Ai > be the first eigenvalue of (4.6) with ( |4.4[ ). We call 
(4.7) T c = Xl/a, 



the critical Taylor number, where a is given by (4.2). 
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As fj, — > 1, equations (4.6) are reduced to the following symmetric linear equa- 
tions: 

dp 
dz 

dp 

(4.8) -^ U ^Wr 

— Aug = Xu r , 

du z du r 



-Au z 
— Au r 



= 0, 

= Xug, 



dz dr 



= 0. 



Let the first eigenvalue Ao > of (4.8) with (4.4) have multiplicity m > 1, the cor- 
responding eigenfunctions be Vi (i = 1, •• • , to), and the corresponding eigenspace 
be 

Eq = span{i>i | 1 < i < m}. 

We remark here that under conditions (2.4) and (4.1), the condition fi — > 1 can 
be equivalently replaced by 

(4.9) n = (2 + *)/(l-|i), 

for some 5 > 0. In this case the parameter a in ( |4.2[ ) is 

a = ( V 2 ,i)/(l - r? 2 ) ~ «J/2. 



When the conditions (4.1) and ( |4.9[ ) hold true, fi — » 1 and ri — > oo. In this 
case the equations (3.1) are replaced by (4.3), and the linearized equations of (4.3) 
reduces to the symmetric linear system (4.8). For the approximate problem (4.8) 
with (3.2), we use Rq to denote the number R defined by ( 3.23[ ) and (3.24): 

1 



Rt) 



Will 



[(G($,Vi),Vi) + G(^i,*),Vi)ff] 



Here ipi is given by (3.14) with (h,ip) satisfying 
(D 2 - a 2 ) 2 h = \ <p, 
(D 2 - a 2 )<p = -X h, 
h = Dh = 0, tp = at 



M, 



and $ is defined by 
(4.10) 



By (3.6) we have 



(A-AoBo)* = G(Vi,Vi), 

B a <P = P(0,<S>g,<S> r ). 

(G($,^ 1 ),Vi)h = ] 
(G(Vi j *),^i)h = -(G(^i,^i),*) h . 



Hence, we infer from (4.10) that 



Ri) 



IVill 
1 



||Vi|| 2 " ((A " AojBo)$,$) ' 
We see that A — Ao-Bo i s symmetric and semi-positive definite, and 
G(V>i,^i)l-Ker04 - A B ), $l_Kcr(A - A B ). 
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Therefore it follows that 



Rn — — i 



1 



IV'i 



;((A - A B ) 1/2 $, (A - X B )^H) H < 0. 



On the other hand, it is known that the number R(fi) in (3.231 is continuous on /j, 
and 

R(fi) — > Rq as jj, — i 1. 
Hence we derive the following conclusion. 



Theorem 4.1. For the Taylor problem (3.1)-(3.3) there exist /x < 1 and < 
rjo < 1 such that for any /io < /x < 1 and r/o < r\ < 1 with /i < rj 2 , the parameter 
R = R{p,n) defined by \3.23 ) is negative, i.e., 

R(/j,r))<0, V/i < M < 1, ?7o<»7<1. 



Consequently, the conclusions in Theorem 3.1 hold true. 



5. Formation of the Taylor vortices and Structural Stability 



Assertions (5) and (6) in Theorem 3.1 provide an asymptotic structure of the 
solutions in the physical space when the gap r-i — t\ is small, as observed in the 
experiments. However, for general parameters r\ and /x we can not give the precisely 
theoretic results, and only present some qualitative description. Here we consider 
two general cases as follows. 



Case: /i > 0. Following [151 H5], f° r the eigenvector (h, ip) of (3.16) the function 
h can be taken as positive and has a unique maximum point in the interval (77, 1). 
Therefore, for the eigenvectors defined by (3.14) and ( 3.15[ ), the vector fields (i/'z, Vv) 
and ("02, Vv) are divergence-free and have the topological structure as shown in 



Figure 3.2 Hence to obtain Assertions (5) and (6) in Theorem 3.1 for any /x > it 
suffices to prove that 

(5.1) h"{r) ^0 at r = r ,l,7?, 



where ro G (rj, 1) is the maximum point of h. We conjecture that the property (5.1 ) 
for the eigenvector (h, ip) of (3.16 1 is valid for all < /1 < 1 and < r\ < 1. 

Case: /1 < 0. In this case, the situation is different 
that the vector field (ip 



,i\> r 



Numerical results show 
in (|3. 14h has k > 2 vortices in the radial direction, 



called the Taylor vortices, which has the topological structure as shown in Figure 
|5.1| see pQ. This type of structure is structurally unstable. However, as discussed 
in [4], under a perturbation either in space H or in 



H = < (u z , u r , ug) G H\ I I ru z dzdr = > , 
there are only finite types of stable structures. In particular, if the vector field 



(ip z ,ip r ) in (3.14) is Z?-regular, i.e., h(r) satisfies (5.1 ), then there is only one class of 
stable structures regardless of the orientation. For example, when u £ H/(T {J H), 
the asymptotic structure of the solution u(t, Uq) of ( 3.1 )-( 3.3 1 is as shown in Figure 



5.2 



and when uq £ H \ T, the asymptotic structure of the solution u(t,Uo) is as 
shown in Figure 5.3 It is clear that the class of structures illustrated by Figure 
5.2| is different from that illustrated by Figure |5.3| The first one has a cross the 



channel traveling flow in the z-direction and the second one does not have such a 
cross the channel flow. 
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Figure 5.1. (tp z ,tp r ) has k 
vortices in r-direction. 
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Figure 5.2. The stable struc- 
ture with a perturbation in 
space H/(T{JH) . 




Figure 5.3. The stable struc- 
ture with a perturbation in 
space H\T. 




6. Proof of Main Theorems 



6.1. Proof of Theorem |3.1[ We shall prove this theorem in the following several 
steps. 



Step 1. We claim that the problem |3.1|-(3.3| bifurcates from (u, A) = (0, Ao) 
to a circle S 1 which consists of stead states. 



It is easy to see that the problem (3.1) with (3.2) is invariant for the transition 
in the z-direction 



u(z, r, i) — > u(z + zq, r, t) 



for^o G 



Therefore, if uq is a steady state solution of (3.1 1 with (3.2), then for any zq S 



the function uq{z + zq. r) is also a steady state solution. We can see that the set 

Z = {u {z + z ,r)\ zoER 1 } 
is homeomorphic to a circle S 1 in Hi for any uq G H\. Hence, the singular points 



of \i.l\ with (3.2) appear as a circle. 



It is known in [161 I15j that there exist singular points bifurcated from (u, A) = 
(0, Ao). Thus this claim is proved. 

Step 2. Reduction to the center manifold. We shall use the construction 



of center manifold functions to derive the reduced equations of (3.5) given by 

dx 



(6.1) 



dt ,iy ' (V>u,V»u)h ' 
^=/MA)y+ (G . (u) 'fi* )g , 
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where ipix and 0ia are the eigenvectors of (3.9) corresponding to /3i(A) near A = Ao 
with 



(6.2) 



lim ip\\ = ipi (-01 as in (3.14)), 
A^A ' ' 



lim Aia = 0i (0i as in (3.15)), 



A->A 



and V'ia and tp*^ are the dual eigenvectors of ip\x and 0ia satisfying 

lim 1A = ipl (V»* as in pTrfb), 

(6.3) A ^ A °- . - ZZ 

lim Via — ^i (01 as in (3.18)). 

A->-An ' ' 



Let \1/ : Eq — ► Eq be the center manifold function of (3.5) at A = Ao, where 

E Q = span{V>i,V>i} 5 

E^ = {u e H\ (u,0i)h = 0, (u,i>* x ) H = 0}. 

Let uq = xipi + yipi £ E - Then it is easy to check that G(uq) € Eq. Hence, by 
the center manifold approximation formula in |12L II 1] , we find that 



(6.4) 



* = 4>{x, y) + o(|.x| 2 + \y\ 2 ) + 0(/MA)(M 2 + M 2 )), 

- L x <t> = G(0i, Vi)x 2 + G(0i, 0i)y 2 + (G(0i, 0i) + G{i>x^))xy. 



On the center manifold, u = uo + $(uo)- Therefore from (6.1 )-(6.4) we obtain the 
reduced equations of (3.5) to the center manifold as follows: 

flT, T 

-=A(A)x + -(G(0,0i) + G(0 1 ,0)>*) ff 
at p 

y, 



(6.5) 



+ ^(G(0,0 1 ) + G(0 1) 0),Vi> 

+ (|x| 3 + H 3 ) + £ i(A)0(| s | 3 + M 3 ), 

'|=/3 1 (A)y+-(G(0,0i) + G^i,0),0i*) ff 
at p 

+ ^(G(0,0i) + 0(^,0),^*)* 
P 

+ (M 3 + M 3 ) + £2 (A)0(M 3 + M 3 ), 



where p= (0i,0i) ff = (-01, -0*)^, and 



lim e,(A) =0 (i = 1,2). 

A->An 



Furthermore, direct calculation shows that 

G(0i,0;i) = p^ o + p^2, 



G(0i,0i) = P0o-PV 2 , 

G(0i,0l) = -P0O + P02. 
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where P : L 2 (M) 3 — > H is the Leray projection, and 



f 0, 



V'o = - < 



V'2 = - < 



tpo 



-(a 2 hD^h + a 2 hDh - -ip 2 ), 
2 r 

-UpD*h + hDip + -hip), 

- sin2az((.D*/i) 2 - hDD^h), 

- cos 2az(a 2 hDh - a 2 hD„h </? 2 ), 

2 r 

- cos2az(hDip — tpD*h -\ — <ph), 
2 r 

°-{{D,h) 2 + hDD*h), 



V'2 



— - cos2az((D.,h) 2 - hDD*h), 

— sin2az(a hDh — a hD*h ip ), 

2 r 

— sm2az(hDip — ipD*h -\ — <ph). 
12 r 



Thus, (6.4 1 is rewritten as 

(6.6) (A - X Q B)4> = {x 2 + y 2 )P^ + (x 2 - y 2 )P^2 + 2xyP^ 2 . 
Let 

(6.7) cf> = -[(x 2 + y 2 )^ + (x 2 - y 2 )<fe + 2xyfa], 

4>2 = ( — l/2siri2azip Zl cos2az(p r ,cos2az(pg), 
<p2 — (l/2cos2az(J3 z ,sm2az(l) r ,sin2az(j}g). 



Then we deduce from (6.6) and (6.7) that 



fz = <Pz, fr = ¥V, fe = <po, 



and (p z ,(p r ,ipe) satisfies 



(DD* - Aa 2 )ip r + Aa 2 X {— - n)tp g = Aa 2 H 2 + 2aDH 1 , 
(DD* - 4a 2 )(p e + X Q Kip r = H s , 

la 
p z = 0, ipg = 0, <p r = Dp r = at r = rj, 1, 
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where H i , Hi , and H3 are given by 

Hi = a({D,h) 2 -hDD,h) 1 



(6.9) 



Ho = a 2 hDh-a 2 hDfa 



1 



^ 2 , 



1 



i^3 = a(hDtp — tpD*h H — </?/i). 
r 



Based on (6.8 1 we find 

(G(&,&)+G(&,&),Vi)ir = for 

(G(^i,^)+G(^,Vi),^i)h - for 
{G{fa,fa) + G{fa,fa),fa) H =Q, 
(G(^i,^,) + G(^ 2 ,Vi),Vi)h = 0. 



t = 0,2, 
t = 0,2, 



Then, putting (6.7) into (6.5), we deduce that 



dx 
~dt 



(6.10) 



dt 



fax - -x(x 2 + y 2 )(G(fa, fa) + G(fa,fa), ^*)h 

--^-^(GC^^O + G^i,^),^)^ 

2 „ ~ - - - 

- -xy 2 (G(fa, fa) + G(fa,fa),fa) H 
P 

+ o(\x\ 3 + \y\ 3 )+e 1 (X)0(\x\ 3 + \y\ 3 ), 

fay - -y(x 2 + y 2 )(G(fa,fa) + G(fa,fa),fa) H 
P 

- -y(x 2 - y 2 )(G(fa, fa) + G(fa,fa),fa) h 
P 

- -yx 2 (G(fa,fa) + G(fa,fa),fa) H 
P 

+ o(\x\ 2 + \y\ 2 )+e 2 (X)0(\x\ 3 + \y\ 3 ). 



Direct computation yields 

(G(fa,fa) + G(fa,fa),fa) H 



= (G(fa,fa)+G(fa,fa),ri) H 

= (G(fa,fa) + G(fa,fa),fa) H 
= (G(fa,fa) + G(fa,fa),fa) H . 



Hence (6.10) can be rewritten as 

dx 

(6.11) f 

dy 

dt 
where 



= fax + Rx(x 2 + y 2 ) + o(|.t| 3 + \y\ 3 ) + £l (A)0(M 3 + |^| 3 ), 
= fay + Ry(x 2 + y 2 ) + o(|.t| 3 + \y\ 3 ) + £ 2 (A)0(|.t| 3 + \y\ 3 ), 



(6.12) 



R= — {G{fa + fa 7 fa)+G{fa,fa+fa),fa) H . 
P 



On the other hand, we infer from ( |6.6| ) and (6.8) that 

{A - A B)4 = Pfa for i = 0, 2, 
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Hence we find 



$ = -(00 + 02), 

(A - A S)$ = G(V>i, Vi) = Ptfto + ^2. 



Thus the number (6.12) is the same as that in (3.23) 



Step 3. Proof of Assertions (l)-(4). When R < 0, (x,y) = is lo- 

Aq. Therefore u = is a locally 



cally asymptotically stable for (6.11) at A 



asymptotically stable singular point of (3.5). By the attractor bifurcation theorem, 



Theorem 6.1 in p. 153 in [5], the problem (3.1 )-(3.3) bifurcates from (w, A) = (0, Ao) 
to an attractor A\ which attracts an open set U \T, and Assertions (1), (3) and 
(4) hold true. 



In addition, the nonlinear terms in (6.11) satisfy the coercive condition in the 



S ll -attractor bifurcation theorem (Theorem 5.10 in [5]) and the conclusion in Step 
1., and Assertion (2) follows. 



Step 4. Attraction in C-norm. It is known that for any initial value u G H 



there is a time to > such that the solution u(t,uo) of (3.1 )-( 3.3 ) is analytic for 



t > to, and uniformly bounded in C r -norm for any r > 1; see Theorem 1 in [7]. 
Hence, by Assertion (4), for any u € U \ T we have 



.13) 



lim min \\u(t, Uo) — Vo\\c r = 0. 

t-¥oo v €A\ 



Step 5. Structure of solutions in A\. By Assertion (3), for any steady 
state solution uq — (u z ,u r ,ug) e .4a, the vector field u = (u z ,u r ) of uq can be 
expressed as 



(6.14) 



u z = 7 cos a(z + zo)D*h(r) + wi(z,r,f3 1 ), 
u r = a7sina(z + zo)h{r) + W2(z,r,fii), 



for some zo £ R 1 , where 

7 = \fc{X)/R\ l '\ w i = o{\p 1 \ 1 ' 2 ) for* = 1,2. 



As in the proof of Theorem 4.1 in [8], we deduce that the vector field (6.14) is 



D-regular for all < A — Ao < e for some e > 0. Moreover, the first order vector 
field in pi 



(6.15) 



(v z , v r ) — (7COsa(z + zo)D^h(r), a7sina(z + zo)h(r)) 



has the topological structure as shown in Figure [372 
Furthermore, it is easy to check that the space 



H = ^ u = (u z ,u ri ug) G H\ / ru z drdz = 

Jm 
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is invariant for the operator L\ + G defined by (3.4 1. To see this, since u is z-periodic 
and u = at r = 1 , 77, we have 



L /•! 



Jr) 



r(u ■ V)u z drdz 



Jo 


f 1 du z 


drdz 


= -L 


L f 1 d{ru 
Jr, dr 


u 7 a 


- r 


f 1 du z 
I r——u z 


drdz 


Jo 


Jn oz 




= 


VueiJ, 





L i-l 



rAu z drdz = 0, Vu € H. 



>0 Jr] 

Thus we see that H is invariant for L\ + G. 



Therefore, for the vector field (6.13) we have 



ru z drdz = 0. 



M 



By The Connection Lemma and the orbit-breaking method in [6], it implies that 



the vector field (6.14) is topologically equivalent to its first order field (6.15) for 
< A - A < e. 

Step 6. Proof of Assertions (5) and (6). For any initial value w € 
£A(rlJ#), we have 



(6.16) 



U = ^2 a kZk + W , 



fe=l 



where wq G H , and for any k = 1, 2, • • • , 

e fe = (e fc (r),0,0), 



and efe(r) satisfies that 



e k (r)dr ^ 0, 



D*De k = -p k e k , 
< pi < 92 < ■ ■ ■ 



Make the decomposition 



E = span{ei,e 2 , • • •}. 
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Then the equation (3.5| can be decomposed into 



.17) 



(1(2 

~dl 
dw 

~dt 



= L x e 



L x w + G{w) 



ee£, 

w e Hi , 



(e(0),io(0)) = C^2a k e k ,w ). 



It is obvious that 

L\e k = -p k e k . 

Hence for the initial value ( 6.16[ ), the solution u(t, u a ) of (6.17) can be expressed as 

u(t, u ) = ^ a k e~ Pkt e k + w(t, u ), 
(6.18) 



f 



e k (r)dr =/= 0. 



By (6.13) we have 



lim \\w(t,uo) - voile- = 0, v e Ax, 



which implies by Step 5 that w(t,Uo) is topologically equivalent to (6.15) for t > 



sufficiently large, i.e., w(t,uo) has the topological structure as shown in Figure 3.2 



By the structural stability theorem, Theorem 2.2.9 and Lemmas 2.3.1 and 2.3.3 



(connection lemmas) in [6], we infer from (6.18) that the vector field in (6.18) 
is topologically equivalent to either the structure as shown in Figure |3.1[ a) or 
the structure as shown in |3.1[ b), dictated by the sign of a k(j in (6.16) with fc — 
min{fc|afe ^ 0}. Thus Assertion (5) is proved. 

Assertion (6) can be derived by the invariance of H under the operator L\ + G 
and the structural stability theorem with perturbation in H, in the same fashion 
as in the proof of Theorem 2.2.9 in |^ by using the Connection Lemma. 

The proof of Theorem |3.1| is complete. 



6.2. Proof of Theorem 3.2 , When R > 0, by Theorem A. 2 in [TO], we infer from 



the reduced equation (6.11) that the transition of (3.1)- (3.3) is of Type-II. In the 



following, we shall use the saddle-node bifurcation theorem, Theorem A. 7 in [TU], 
to prove this theorem. Let 

H* = {(u Zl u r ,u e ) € H\ u z (-Z,r) = -u z {z,r)}, 

H* =H 1 f]H*. 

It is easy to see that the space H* is invariant under the action of the operator 
L x + G defined by Q: 

(6.19) L X + G:H{^H\ 

and the first eigenvalue /?i(A) of Lx : HI — > H* at A = \q(T = T c ) is simple, with 



the first eigenvector ipi given by (3.14). Hence, the number R in (6.12) is valid for 



the mapping (6.18), i.e., 

(G(xV>i + *(*), VD = Rx 3 + °(M 3 )- 

Thus, it is readily to check that all conditions in Theorem A. 7 in [TU] are fulfilled 
by the operator (6.19). By Step 1 in the proof of Theorem 3.1 each singular point 
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of (6.19) generates a singularity circle for L\ + G in H. Therefore, Theorem 3.2 
follows from Theorem A. 7 in [ID] . 

The proof of Theorem |3.2| is complete. 
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